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Abstract 

We study the light CP-even neutral Higgs boson production in association with an 
electron and a jet at the possible CERN large hadron-electron collider within the minimal 
supersymmetric standard model. We investigate the possible supersymmetric effects on 
this process and compare our standard model numerical results with those in previous 
work. We present the leading-order and QCD next-to-leading-order corrected total cross 
sections and the distributions of the transverse momenta of the final electron, the light 
neutral Higgs boson, and jet in the minimal supersymmetric standard model. Our results 
show that the scale dependence of the leading-order cross section is obviously reduced by 
the QCD next-to-leading-order corrections. The K factor of the QCD correction to the 
total cross section at the large hadron-electron collider varies from 0.893 to 1.048 when 
the factorization/renormalization scale /i goes up from O.lmz to 3.8mz in our chosen 
parameter space. 

PACS: 12.60.Jv, 12.38.Bx, 14.80.Cp 
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I. Introduction 



One of the most significant tasks for high-energy experiments is to search for scalar Higgs 
particles [THS]. Although the standard model (SM) [6j has achieved impressive experimental 
success, the Higgs boson, which is predicted by the SM for spontaneous electroweak symmetry 
breaking, remains a mystery. Moreover, there exists the problem of the quadratically divergent 
contributions to the corrections to the Higgs boson mass, which is the so-called naturalness 
problem. Alternative conceptional difficulties, such as the hierarchy problem, the necessity of 
the tuning and the nonoccurrence of gauge coupling unification at high energies, suggest that 
the SM is probably the low-energy limit of a more fundamental theory. 

As the most hopeful extensions of the SM, the supersymmetric (SUSY) models can solve 
such problems mentioned above. The minimal supersymmetric standard model (MSSM) |71IH] is 
the simplest one among all the SUSY extensions of the SM. In this model, two Higgs doublets Hi 
and H2 give masses to up- and down-type fermions. The Higgs sector consists of three neutral 
Higgs bosons, one CP-odd particle two CP-even particles and and a pair of 

charged Higgs bosons [H^). However, these Higgs bosons have not been directly explored 
experimentally until now. The LEP experiments provided lower mass bounds as: for the SM 
Higgs boson mno > 114.4 GeV (at 95% C.L.), and for the MSSM bosons rriho > 92.8 GeV and 
niAo > 93.4 GeV for tan/3 > 0.4 (at 95% C.L.) 

Recently, a possible high-energy collider in e~-p collision mode at the LHC, the large hadron- 
electron collider (LHeC), has been sketched [TT|IT2]. There will exist a rich physics program |13j . 
The LHeC can be used to accurately determine the parton dynamics and the momentum 
distributions of quarks and gluons in the proton, and furthermore it may play a significant role 
in the discovery and interpretation of new physics. The incoming proton beam at the LHeC has 
an energy Ep = 7 TeV and the energy of incoming electron is considered as E^, = 50 — 200 GeV 
according to several scenarios, with the center-of-mass system energy of -y/i = 2 a/ EpEe ~ 
1.18 — 2.37 TeV. It seems that the LHeC provides a cleaner environment than a hadron-hadron 
collider in accessing the couplings of the Higgs boson to gauge bosons. 

The production channel e~p — )■ e'hPj + X , a. neutral current (NC) process at the LHeC, 
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attracted the physicist's attentions. In Ref. [H] it is pointed out that the electron reconstruction 
in the NC process is superior with respect to that of the missing neutrino in the charged current 
process, e~p — )■ Ueh'^j + X, and the NC process has the potential to increase the overall Higgs 
boson signal efficiency, and there they studied the use of forward jet tagging as a means to 
secure the observation of the Higgs boson in the — ?■ bb decay mode and to significantly 
improve the purity of the signal. The QCD next-to-leading-order (NLO) corrections to the SM 
Higgs productions of e~p — > e^H^j + X and e~p — )■ VeH^j + X processes at the LHeC were 
calculated by B. Jager in Ref. [12]. Moreover, not only does this channel provide a spectacular 
signature {e~bbj), but also the lightest Higgs production in MSSM via vector boson fusion 
with unusual visible decays is possible [TB]. The coupling strength of the lightest Higgs 
with Z^Z^ is distinguished from the SM Higgs one with an additional factor sin(/3 — a), where 
P is related to the ratio of the vacuum expectation values and a is the mixing angle of the 
two CP-even Higgs states. Therefore, we may disentangle between the SM Higgs and the light 
MSSM CP-even Higgs by measuring the cross section for e~p — )■ e~h^j + X at the LHeC when 
I sin(/3 — a)| is smaller than 1. Besides, in order to find new physics it requires sufficiently precise 
predictions for the new physics signals and their backgrounds with multiple final particles which 
cannot be separated in experimental data entirely. Therefore, the higher order QCD predictions 
for these reactions are necessary. 

In this paper, we calculate the full QCD NLO corrections to the process e~p — )■ e~h^j + X 
at the LHeC and estimate the capability of the LHeC to access the light MSSM CP-even Higgs 
boson in the e~h^j production. The numerical results at the leading-order (LO) are compared 
with those in Ref. [2]. The paper is organized as follows: We describe the technical details of 
the related LO and QCD NLO calculations in both the SM and the MSSM in Sees. II and HI, 
respectively. In Sec. IV we give some numerical results and discussions about the QCD NLO 
corrections in the MSSM. Finally, a short summary is given. 
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Figure 1: The LO Feynman diagram for the partonic processes e g — ?• e h'^q {q = 
u, u, d, d, c, c, s, s). 



II. LO cross sections 

The LO and QCD NLO calculations are carried out in 't Hooft-Feynman gauge. The FeynArts 
3.4 package [T7] is adopted for generating Feynman diagrams and subsequently converting 
them to corresponding amplitudes. The FormCalc 5.4 program [18] is applied to reduce the 
amplitudes. 

In calculating the e~p — )■ e^h^j + X process in the MSSM, we neglect the u-, d-, c-, s-quark 
masses {rriu = = rric = rris = 0), and do not consider the partonic processes with incoming 
(anti)bottom-quark due to the heavy (anti)bottom-quark suppression in parton distribution 
functions (PDFs) of proton. That means we involve the contributions of the following partonic 
processes in our LO calculations: 

e"(pi) + g(p2) -> e"(p3) + /i°(p4) + g(p5), (g = m, m, rf, J, c, c, s, s), (2.1) 

where Pi{i = 1,...,5) represent the four-momenta of the incoming electron, partons, and the 
outgoing electron, /i°-boson and jet, respectively. The LO Feynman diagram for the partonic 
processes (12.11) is depicted in FigHJ 

The expression of the LO cross section for the partonic process e~g — )■ e~h^q can be written 
in the form as 

aLois,e~q ^ e~h°q) = ^ j ^\MLo\^dn3, {q = u,u,d,d,c,c,s,s), (2.2) 

where s is the partonic center-of-mass (cm.) energy squared, the summation is taken over 
the spins and colors of final states, and the bar over the summation means taking average 
over the intrinsic degrees of freedom of initial particles, and dQ^ is the three-body phase-space 
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element for the e~g — )■ e~h^q subprocess. Mlo in Eq. f|2.2p is the tree-level amplitude for the 
partonic process e^q — > e~h^q . The coupling of the SM Higgs boson (H^) to the pair can 
be expressed as g^zz ~ ^l^fi'^^? while the light CP-even SUSY Higgs boson to the pair 
is expressed as g^^zf^^ = sin(/3 - a)g^"'. 

The LO total cross section for the e^p — )■ e^h^j + X process at the LHeC can be expressed 

as 



There /x/ is the factorization scale, s is the total cm. energy squared of the electron-proton 
collision, x describes the four-momentum fraction of parton q in an incoming proton with 
the definitions of x = ^, and P is the four- momentum of the incoming proton. Gqjp {q = 
u, u, d, d, c, c, s, s) represent the PDFs of parton q in proton p. 



III..1 Virtual corrections 

In order to compare the results in the MSSM with those in the SM we present the QCD NLO 
calculations in both models. In the NLO calculations, we adopt the dimensional regulariza- 
tion in D = 4 — 2e dimensions to isolate the ultraviolet (UV) and infrared (IR) singularities. 
The wave functions of the external fields are renormalized under the on-shell renormalization 
scheme. The virtual correction to the subprocess e~g — t- e~h^q involves both soft and collinear 
IR singularities. In our calculation we introduce the following counterterms for related wave 
functions in the SM and the MSSM: 



The wave-function renormalization constants of the massless quarks (g = u, d, c, s) in the SM 
are written as 




III. QCD NLO corrections in the MSSM 




(g = c, s). 



(3.1) 




3^ 



(3.2) 
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Figure 3: The representative pure SUSY QCD one-loop Feynman diagrams for the partonic 
process e~q — )■ e~h^q, where q = u,d,c, s and the lower indexes k = 1,2. 

where Af/y = l/e[/y— 7£;+ln(47r) and A/r = l/e//?— 7£;+ln(47r). The explicit expressions for the 
one-loop QCD wave-function renormalization constants of the massless quarks (g = u, d, c, s) 
in the MSSM have the forms as 

^^MSSM ^ _^ [^^^ _ ^ 2a. ^^^^^^ ^_ ^ ^^^^^ 

5^M55M ^ _^ [^^^ _ + ^ [^^^0^ ^2^ gi^2 ^_ ^ ^^^0^ ^2^ ^2g ^^^2 ^3 4) 

where the definitions for the two-point integrals are adopted from Ref. [19], and 9q is the mixing 
angle of scalar quarks {qi, qn), 

qi = qiCosOg - q2sin6q, q^ = qi sin 6q + q2 cos 6q. (3.5) 

The one-loop level Feynman diagrams include self-energy, vertex, box (4-point) and counterterm 
Feynman graphs. We depict the SM QCD vertex diagram in FigJJl and the representative pure 
SUSY QCD (pSQCD) one-loop diagrams are drawn in Figj3l 
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III. .2 Real gluon and light-(anti)quark emission corrections 

The relevant real emission partonic processes can be grouped as (1) e~q — )■ e~h^qg, (2) e~g — >• 
e~h^qq. There the quark notation, g, represents m— , m— , d—, d—, c— , c— , s— and s— quark. The 
real gluon/light-(anti)quark emission partonic channels (1) and (2) at the tree-level contain 
soft and collinear IR singularities. After the summation of the virtual corrections with all the 
real parton emission corrections, the numerical result is soft IR-safe, while collinear divergences 
still remain. It will be totally IR-safe when we include the contributions from the collinear 
counterterms of the PDFs. The IR finiteness can be verified numerically in our numerical 
calculations. 

The IR singularities of the real parton emission subprocesses can be isolated by adopting 
the two cutoff phase-space slicing method [20] . In FigsH] and [5] we present the Feynman di- 
agrams for the real gluon emission subprocess e~(pi)g(p2) — ^ ^ {pz)h^ {P'i)l{p^)9{pQ) and real 
light- (anti) quark emission subprocess e~ {pi)g{p2) — ^ e~(p3)/i°(p4)g(p5)g(p6), respectively. In 
adopting the two cutoff phase-space slicing method we introduce an arbitrary small soft cutoff 
5s to separate the 2 — )■ 4 phase-space into two regions, Eq < 5s\/~^/2 (soft gluon region) and 

> 6sy/^/2 (hard gluon region), and another cutoff 6c to decompose the hard region into 
a hard collinear (HC) region with P2{Pb)-P6 < 6cs/2 and hard noncollinear (HC) region with 

Then the cross sections for the real emission subprocesses e~{q,g) — )■ e^h^q{g,q) can be 
written as 

= + = + + a^. (3.6) 

IV. Numerical Results and Discussion 

In our numerical calculations we take one-loop and two-loop running as in the LO and NLO 
calculations, respectively [S]. The QCD parameters are taken as Nf = 5, A^^ = 165 MeV 
and A*^'^ = 226 MeV. We take the renormalization and factorization scales to be a common 
value as /i = /ir = /i/ and choose the energy scale to be at the mass (i.e., /x = /iq = rriz) by 
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Figure 4: The tree-level Feynman diagrams for the gluon emission partonic process e g — ?• 
e~h^qg {q = u, u, d, d, c, c, s, s). 
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Figure 5: The tree-level Feynman diagrams for the light-quark emission partonic process e~g — )■ 
e~hPqq (g = m, c?, c, s). 

default, which characterizes the typical momentum transfer in the process e~p — ?■ e^h^j + X . 
The relevant SM parameters are taken as rrie = 0.511 MeV, rrih = 4.2 GeV, rrit = 171.2 GeV, 
mw = 80.398 GeV, mz = 91.1876 GeV and Gp = 1.16637 x 10"^ GeV""^ [9\, and thus we get 
a = 1/132.34 by adopting the relation of a = ^G^m^s^. We use the PDFs of CTEQ6L1 
and the CTEQ6M in the LO and NLO calculations, respectively [21j. 

The related SUSY parameters, such as the mixing angle of the MSSM Higgs fields a and 
masses of the light CP-even neutral Higgs boson, gluino, and scalar quarks, are obtained from 
the FormCalc package, except otherwise stated. The input parameters for the FormCalc 
program are Ms, M2, Af, m^o, /i and tan/3. There Mq = Mu = Mo = Ms and the soft 
trilinear couplings for squarks q being equal, i.e., Aq = Ai = Aj are assumed, and the grand 
unification theory relation Mi = (5/3) tan^ 9wM2 is adopted for simplification. In our numerical 
calculation, we set Ms = 400 GeV, M2 = 110 GeV, m^o = 150 GeV, fi = -200 GeV, 
Af = 800 GeV, tan/3 = 3, and nig = 230 GeV in defauh. Then we get sin(/3 - a) = 0.9347, 
rui^ = 198.17 GeV, m^^ = 579.67 GeV, ruu, = rUc, = 397.07 GeV, ruu^ = m^^ = 398.76 GeV, 

= e-c = 7r/2, mj^ = m-s^ = 400.62 GeV, mg^ = m,-^ = 403.52 GeV, 9^ = 9^ = and 
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TUho = 98.36 GeV. In the FormCalc program the radiative corrections to the MSSM Higgs 
boson masses up to two-loop contributions are involved, and the expressions related to the 
mass of the light neutral CP-even Higgs boson in Ref. [22] are adopted, where the input 
parameters of m;,, m^, m^^ and m^^ are necessary. 

The verifications for the total QCD NLO correction being independent of the two cutoffs Sg 
and 6c are made. We calculate the total QCD NLO corrections to the e~p — > e^h^j + X process 
in the MSSM at the LHeC with the cutoffs Ss running from 10"^ to 10"^ 4 = 5^/200, and 
jj, = jjLQ = VTiz- The results show that although the three-body correction [Aa*^^-' = cr^+cr'^+o"'^'"] 
and four-body correction [Acr*^^^ = a^'~^\ depend strongly on the cutoff (5c), the final total 
QCD NLO correction /\aNLOi which is the summation of the three-body and four-body terms, 
i.e., Aaj\iLo = Aa*^^^ + Acr*^^) is independent of the two cutoffs within the statistic errors. The 
independence of the full QCD NLO corrections to the e~p — )• e'q{q = u^u) ^ e~h^j+X process 
on the cutoffs 6s and 6c provides an indirect check for the correctness of the calculations. In 
further numerical calculations, we fix 5^ = 8 x 10~^ and 6c = 5^/200. 

We made the comparison of our LO numerical results for the process e~p e~H^j + X in 
the SM at the LHeC with the corresponding results read out from Fig.2 in Ref. [H], and find 
that they are coincident with each other within the statistic errors. 

In the following LO and NLO numerical calculations, we adopt the massless four-flavor 
scheme and put the restriction of > p^j on the jet transverse momentum for one-jet 
events. For the two-jet events (originating from the real corrections), we apply the jet al- 
gorithm in the definition of the tagged hard jet with R = 1, i.e., if final state two partons 
satisfy a/ Ar^^ + A^^ < 1 (where Arj and A0 are the differences of rapidity and azimuthal angle 
between the two jets), we merge them into a single jet. We use the so-called "inclusive" scheme 
and keep events with one or two jets. We require that there is one jet with > p^j, and set 
p^^j = 30 GeV by default in following calculations. Furthermore, to reduce the background of 
the Higgs signals, we require the final electron with the following cuts 

> 30 GeV, \7]^\ < 5. (4.1) 

We plot the dependence of the LO and QCD NLO corrected total cross sections for the 
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e~p — )■ e~h^j + X process in the MSSM on the renormahzation/factorization scale /i in FigJHI^a). 
The corresponding K factor defined as = versus the energy scale is presented in 

FigJEl^b). Figur^a) shows that the LO curve is obviously dependent on the energy scale /z, 
although there only the factorization scale is involved in the convolution with the PDFs of the 
initial parton. If we define the scale uncertainty parameter as = |^~^^^|^| in the scale 
range of [/ii = |yU0)/^2 = 3/io], from FigJHI^a) we can get the uncertainty parameters for the LO 
and QCD NLO corrected cross sections having the values as rjio = 10.35% and rj^Lo = 1-58%, 
respectively. It is obvious that the dependence of the LO total cross section on the energy scale 
is significantly reduced by the QCD NLO corrections. We can read out from FigJH^b) that 
when the energy scale varies from 0.2m^ to S.Sm^, the value of the K factor increases from 
0.893 to 1.048. In the following, we choose fi = fiQ except otherwise stated. 




Figure 6: (a) The dependence of the LO and QCD NLO corrected total cross sections for the 
e~p — )■ e~/i°j + X process on the renormahzation/factorization scale /i = /x^ = /^/ in the MSSM, 
where we take Ep = 7 TeV and = 140 GeV. (b) The corresponding K factor of Figl6]^a) 
versus the energy scale n (where we define K = ^^)- 

We plot the LO and QCD NLO corrected total cross sections for the e^p e'h^j + X 
process in the MSSM as a function of the incoming electron beam energy E^. running from 
50 GeV to 200 GeV in FigJTl^a), that corresponds to the cm. colliding energy range of ^/s ~ 
1.18 — 2.37 TeV. The corresponding K factors are depicted as a function of the incoming 
electron beam energy E^. in FigJT^b). In FigJTI^a) the full line is for the QCD NLO corrected 
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total cross section for the e~p — )■ e~h^j + X process, and the dotted hne for the LO cross 
section. We can see from FigsJTl^a) and[71^b) that the QCD NLO corrections reduce shghtly the 
LO total cross sections for the process e~p e^h^j +X in the plotted incoming electron beam 
energy range, and the production rate increases with E^.. In FigjTl^c) we depict the K factor 
versus electron beam energy Ee, the energy scales /i being O.SyUo and 3yUo separately. We can 
see from Fig|71^c) that the K-factor uncertainty, A/T = = 3fio) — K{fj, = 0.5/io), ranges 
from 12.79% to 7.13% when E^ goes up from 50 GeV to 200 GeV. 

I 1 ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 1.03 I 1 . 1 . 1 . 1 . 1 . 1 . 1 1 1 




60 80 100 120 140 160 180 200 60 80 100 120 140 160 180 200 



Ee(GeV) Ee(GeV) 




Figure 7: (a) The dependence of the LO and QCD NLO corrected total cross sections for 
the e~p — 7- e^h^j + X process in the MSSM on the incoming electron beam energy Ee in the 
MSSM, where we take fi = fiQ, Ep = 7 TeV, and E^ = 140 GeV. (b) The corresponding K 
factor [K = ) versus the incoming electron beam energy Ef.. (c) The K factor versus the 
incoming electron beam energy E^ with /x = 0.5/io and /i = 3/io, respectively. 

The curves for the LO and QCD NLO corrected cross sections for the process e~p — >• 
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e~h^j + X as a function of tan /3 are drawn in FigJSI^a), where the corresponding values of 
m/jO are also shown on the x axis in FigsJHl^a) and El^b). The values of m^o and of the other 
parameters are those given above. In FigIH](a), we can see that both curves go down rapidly in 
the region of 2 < tan/3 < 6 (85.52 GeV < rriho = 113.14 GeV). Then the curves go up slowly 
after the values reach their corresponding minimal values at position around tan/3 ~ 7.5. The 
relevant K-factor {K = (Tnlo/o'lo) versus tan/3 (and m/jo) is plotted in FigJHI^b). The K 
factor generally has a constant value of about 0.99. We further depict two curves for the K 
factors with fi = 0.5/io and fi = 3fiQ separately, as a function of tan/3 (and 171^0) in FigJHI^c). 
We can read out from FigJHI^c) that the K-factor uncertainty due to the scale n, defined as 
AK = K{fi = 3yUo) — = 0.5yUo)) is in the range from 4.06% to 6.29% when tan/3 {rjif^o) 
varies from 2 (92.80 GeV) to 50 (121.64 GeV). 

For the comparison of the results for the processes e~p — t- e~H^j + X in the SM and 
e-p e'h^i+X in the MSSM at the LHeC, we read out the data in the MSSM from Fig|S][a) 
at the positions of tan/3 = 3,7,18,38 respectively, and list these results together with the 
corresponding SM ones in Table [U All SM parameters, including the mass of the SM Higgs 
boson, have the same values in both the SM and the MSSM calculations. The relative difference 
between the cross sections in both models, is defined as 8^10 = sm x 100%. These 

"nlo 

numbers are obtained by adopting the values of the renormalization/factorization scale and 
the input SUSY parameters mentioned above. From the table we can see that the relative 
difference, Snlo, between the cross sections in both models can reach the value of —10.57% up 
to the QCD NLO, when we take tan /3 = 3 for the MSSM. 

In Figjni^a) we depict the LO and QCD NLO corrected cross sections for the process e~p — >■ 
e'h^j + X as a function of mass m^o (and nif^o). As we know, the light CP-even Higgs boson 
mass m/jO depends on the CP-odd Higgs boson mass m^o, when the other related SUSY and 
SM input parameters are fixed. The values of m^o corresponding to different m^o values are 
also shown on the x axis in Figsin](a) and [H^b). In FigJHI^a), we see that the cross sections 
increase rapidly in the range of 100 GeV < rriAO < 180 GeV (It corresponds to the range of 
80.92 GeV < rrif^Q < 106.02 GeV). After reaching their maximal values at position of m^o = 
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Figure 8: (a) The LO and QCD NLO corrected total cross sections for the e~p — )■ e~hPj + X 
process as a function of tan^ {tuao fixed) and the corresponding mass of the hght CP-even 
neutral Higgs boson ruho in the MSSM, where we take /x = /xq, Ep = 7 TeV, and = 140 GeV. 
(b) The corresponding K factor {K = versus tan/3 and m^o. (c) The K factor versus 

tan/3 and m/jO with /x = 0.5/xo and /x = 3/xo, respectively. 
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Table 1: The numerical results of the a^^^^^^ , (^nlo^^ for tan/3 = 3, 7, 18, 38 obtained 
from FigJHl^a), and the corresponding SM results of the crfoNLO process 
e^p — i- e^H^j + X are listed in the table, where we take the same SM parameters 
and the mass of the Higgs boson (m/jo = niHo) in both the SM and the MSSM 

„MSSM _„SM 

calculations. 5nlo is defined as '"^^o ^.^/^vlo x 100%. 

'^NLO 



220 GeV, the LO and QCD NLO corrected cross sections decrease gently. The corresponding 
K factor versus m^o (and m/jo) is displayed in Fig|9]^b). The K factor seems to be stable and 
has the value around 0.99. We can see that when we fix the energy scale /i = /iq, the QCD 
NLO corrections in the MSSM generally reduce the LO cross section by about 1%, while the 
pure SUSY QCD (pSQCD) NLO contributions are negligibly small and the relative pSQCD 
corrections have the values less than 0.01%. But as seen earlier [Fig|6]^b)], when fi = 0.2m^ 
(3.8mz) the QCD NLO relative correction reaches -10.7% (4.8%). 




Figure 9: (a) The LO and QCD NLO corrected total cross sections for the e~p — )■ e~h^j + X 
process as a function of the masses of the CP-odd Higgs boson and the light CP-even 
neutral Higgs boson {tan (3 fixed) in the MSSM, where we take /x = /io, Ep = 7 TeV and 
Ee = 140 GeV. (b) The corresponding K factor(i^' = (Jnlo/(^lo) as a function of m^o and 
rriho. 
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The distributions of the transverse momenta of the final particles at the LO and up to the 
QCD NLO, and their corresponding K factors for the process e~p — )■ e~h^j + X are depicted in 
FigsIIOi;a,b,c), where we define K = ^^^/^^- In FigsHDJ^a), (b) and (c), the distributions 
of transverse momenta and K factors are for the final electron, the light CP-even neutral Higgs 
boson and jet, respectively. We can find that there is no obvious distortion induced by the 
QCD NLO corrections for the and p^" distributions, while the shape distortion for the p!^^ 
distribution is not negligible since the K factor of the p!^^ distribution varies in the range of 
0.865 < K..t < 1.049. 

V. Summary 

In this paper we calculate the full QCD NLO corrections to the light CP-even neutral Higgs 
boson production associated with an electron and a jet in the MSSM at the possible CERN 
LHeC. We investigate the uncertainty of the integrated cross sections induced by the factoriza- 
tion/renormalization scale, and present the LO and QCD NLO corrected total cross sections 
and the distributions of the transverse momenta of final particles. By adopting the definition of 
the scale uncertainty parameter in the scale range of [/ii = |/io,Ai2 = 3/io] as 77 = I ^^^^^^^^y^ I > 
we obtain the scale uncertainty parameters for the LO and NLO corrected cross sections are 
10.35% and 1.58%, respectively. It is clear that the scale dependence of the LO cross section is 
obviously improved by the QCD NLO corrections. We find that there is no obvious distortion 
induced by the QCD NLO corrections for the p'j- and p^" distributions, and the K factor of 
the QCD correction to the total cross section at the LHeC varies from 0.893 to 1.048 when the 
factorization/renormalization scale /x goes up from 0.2m^ to 3.8m^ in our chosen parameter 
space. 
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Figure 10: (a) The LO and QCD NLO corrected differential cross sections and the corre- 
sponding K factor {k = ^^^/^^ for the process e'p e'h^j + X . (b) The LO and QCD 

NLO corrected differential cross sections ^ and the corresponding K factor (^K — / ^^^^ 
for the process e~p — >■ e~hPj -\- X . (c) The LO and QCD NLO corrected differential cross sec- 
tions -0^ and the corresponding K factor {^K = ''^^^f /^fg-^ for the process e~p e^h°j + X 
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